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2 Abstract 

3 A revised version of the Microphysics of clouds with Relaxed Arakawa-Schubert and 

4 Aerosol-Cloud interaction (McRAS-AC), including, among others, the Barahona and Nenes 

5 ice nucleation parameterization, is implemented in the GEOS-5 AGCM. Various fields from a 

6 10-year long integration of the AGCM with McRAS-AC were compared with their 

7 counterparts from an integration of the baseline GEOS-5 AGCM, and with satellite data as 

8 observations. Generally using McRAS-AC reduced biases in cloud fields and cloud radiative 

9 effects are much better over most of the regions of the Earth. Two weaknesses are identified 

10 in the McRAS-AC runs, namely, too few cloud particles around 40S-60S, and too high cloud 

1 1 water path during northern hemisphere summer over the Gulf Stream and North Pacific. 

12 Sensitivity analyses showed that these biases potentially originated from biases in the aerosol 

13 input. The first bias is largely eliminated in a sensitivity test using 50% smaller aerosol 

14 particles, while the second bias is much reduced when interactive aerosol chemistry was 

15 turned on. The main drawback of McRAS-AC is dearth of low-level marine stratus clouds, 

16 probably due to lack of dry-convection, not yet implemented into the cloud scheme. Despite 

17 these biases, McRAS-AC does simulate realistic clouds and their optical properties that can 

1 8 improve with better aerosol-input and thereby has the potential to be a valuable tool for 

19 climate modeling research because of its aerosol indirect effect simulation capabilities 

20 involving prediction of cloud particle number concentration and effective particle size for 

21 both convective and stratiform clouds is quite realistic. 

22 

23 1 Introduction 

24 Traditionally, meteorologists focused on severe weather and precipitation forecasts. Not 

25 much attention was paid to cloud water. There are two reasons for this. First, in-cloud water is 

26 generally less than 5% of precipitation generated in a typical weather episode; second, 

27 weather forecast is useful for a week or less, and on that time-span, cloudiness and its 

28 radiative effects on the synoptic weather systems are small. Consequently, ad hoc ways to 

29 assess cloud radiative forcing were deemed sufficient. However, once the emphasis of 

30 forecasting turned to climate, radiative forcing and everything that affects it including 

3 1 influence of aerosols on clouds, cloud radiative effects (CRE), and greenhouse gases become 

32 very important. Among them, cloud-aerosol interaction (Andreae and Rosenfeld, 2008) is in 
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1 early stages of development (e.g., Quaas et al., 2003, Roelofs et al., 2006, Sud and Lee, 2007; 

2 Morrison, and Gettelman 2008, Liu et al., 2011). Our paper deals with evaluating a recently 

3 revised and updated aerosol-cloud interaction (AC) module of an established cloud scheme 

4 called McRAS (Microphysics of clouds with Relaxed Arakawa-Schubert convection) that is 

5 implemented in the GEOS-5 AGCM for evaluating its performance. 

6 The pioneering works of Gibbs (1876, 1878), and Kohler (1936) laid the foundation of the 

7 physics of aerosol-activation from first principles of thermodynamics. Recent developments 

8 of physically based aerosol activation parameterizations are due to Abdul-Razzak and Ghan 

9 (2000, 2002), Nenes et al. (2001), Nenes and Seinfield, (2003), Liu and Penner (2005) and 

10 Barahona and Nenes (2009). More aerosols generally increase the number density of cloud 

11 particles (CP) (Twomey, 1979; Seinfeld and Pandis, 1996) and thereby suppress 

12 autoconversion and accretion that fonn precipitating hydrometeors (e.g., Albercht, 1989; 

13 Seifert and Beheng, 2001, 2005). We cite a few recent studies that show the impact of 

14 aerosols on i) weather and climate prediction (Krishnamurti et al, 2009; Sud et al, 2009; 

15 Wilcox et al, 2009) with GEOS-4 AGCM; ii) the diurnal and seasonal cycles of precipitation 

16 (Kim et al., 2010); iii) the weekly cycle of precipitation over central North America (Bell et 

17 al., 2009b); iv) increase in the incidence of tornados (e.g., Rosenfeld and Bell, 2011) and 

18 lightening (Bell et al, 2009a), v) the vertcal stability of the atmosphere leading to an elevated 

19 heat pump hypothesis affecting the Indian monsoons (Lau and Kim, 2007), and vi) freezing of 

20 in-cloud drops with release of latent heat of freezing (Rosenfeld et al., 2000, 2006). The 

21 debate on how and/or how much do aerosols influence different clouds to foster or suppress 

22 precipitation yield from different cloud types continues (e.g., Koren et al., 2012; Gunturu, 

23 2010, and Li et al., 2011). An outcome of increasing cloud particle number concentrations 

24 (CPNC) is the accompanying reduction in cloud particle sizes in the distribution. It slows the 

25 auto-conversion and accretion of cloud drops to form precipitation size hydrometeors and that 

26 gives liquid cloud particles and embryonic raindrops the time to ascend in the convective 

27 updrafts and glaciate at subfreezing temperatures to latent heat of freezing and further boost 

28 the updraft buoyancy. How this plays out in the real world under a variety of convective 

29 scenarios, depends upon how much further the convective towers ascend and how much 

30 additional condensate and precipitation is generated. 

31 Numerical models can simulate all of the above features, if the cloud-physics processes 

32 are realistically parameterized and the amdient atmosphere has realistic aerosols. However, 
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1 many GCMs still obtain clouds invoking a number of simplifying and/or ad hoc assumptions 

2 (Bougeault and Geleyn, 1996; Randall 2010) that often ignore the aerosol effects on clouds. 

3 At the present time, even understanding of the climatic impact of aerosols remains uncertain 

4 (IPCC 2007), while modelers endeavor to include them in their cloud schemes and 

5 intercompare their performance with other models (Bellouin et ah, 2011). Without realistic 

6 aerosol input, simulated CPNC consisting of liquid and ice cloud particle number 

7 concentrations (hereafter LPNC, IPNC), also get affected and real benefits of including 

8 aerosol-cloud interaction become uncertain. A notable inference of cloud seeding 

9 experiments, as summarized by Cotton and Pielke (1995), is that with a limited window of 

10 opportunity, cloud seeding becomes a hit or miss venture and its full potential is not realized. 

1 1 One easily infers that to observationally verify the impact, the atmospheric aerosols too must 

12 be known reasonably well at the scales of clouds. 

13 Aside from CPNC enhancement via shattering and splintering of the already existing 

14 cloud particles aerosol activation for liquid particles or ice nucleation for cloud ice partcles is 

15 the only source of new CPs. A well designed aerosol-cloud-radiation interaction model, such 

16 as McRAS, can therefore provide a better physical basis for detennining the benefits of 

17 interactive aerosols in a cloud scheme. McRAS was renamed as McRAS-AC after the aerosol 

18 cloud interaction was included. The question we want to address is whether McRAS-AC is a 

19 worthwhile option to simulate the cloud optical properties and climate in the GEOS-5 GCM? 

20 Some issues with simulating mixed phase clouds are proverbial. The major one is getting 

2 1 mixed phase water and ice mass fractions, effective sizes, and apportionment of precipitation 

22 in to liquid and ice partcles. Based on the vertical velocity and entrainment rate, IN nucleate 

23 cloud ice particles. In addition we perform mass transfer from liquid to ice particles by 

24 Bergeron-Findeisen process in the following manner. At subfreezing temperatures, the vapor 

25 pressure differences over cloud water and ice particles are large enough to produce a 

26 substantial vapor pressure gradient between them to induce mass transfer of cloud water-to- 

27 cloud ice through the intervening atmosphere. Simultaneously, precipitating hydrometeors 

28 collect cloud water/ice particles during fall through a cloud (even from clear air that happens 

29 to be supersaturated with respect to the ice). In this way, precipitation removes cloud particle- 

30 mass and reduces the in-cloud CPNC. Altogether, these processes add up to make the sink 

3 1 term of CPs for which each of the components must be parameterized. The sum of source and 

32 sink terms yields the time rate of change of mass and number concentration as shown in 
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1 Morrison and Gettelman (2008). To close the system, one also needs a precipitation 

2 microphysics scheme. Thus an end-to-end aerosol-cloud interaction parameterization starts 

3 with aerosols activating as CCN and/or IN nucleating ice to receive condensation and/or ice 

4 deposition; needs a reasonable treatment of mass transfers among liquid, ice and vapor phases 

5 of cloud water with precipitation microphysics for liquid and ice clouds. Most present-day 

6 GCM modelers have started to include prognostic parameterizations of the direct and indirect 

7 aerosol effects in the cloud schemes, but several complexities and uncertainties still confound 

8 them. For example, aerosol input, implying their size distribution and speciation and their 

9 hygroscopic properties, and activation/nucleation potential for CCN/IN, which must be 

10 inferred from aerosol chemistry, are themselves uncertain and are going through extensive 

1 1 validation and upgrading tests. The ultimate challenge is to make the aerosol-cloud-radiation 

12 scheme realistic enough to affect the key aspects of cloud -radiative interactions realistic 

13 enough to better simulate the climate, and respond to aerosol direct and indirect effects. 

14 Section 2 gives a brief description of GEOS-5 AGCM hosting McRAS-AC as one of its 

15 options. Section 3 gives the simulation experiments. Section 4 gives simulation results and 

16 key biases. Section 5 has summary and conclusions and research directions to make MCRAS- 

17 AC simulations more realistic. 

18 

19 2 Cloud Schemes: GEOS-5 GCM and McRAS-AC 

20 2.1 GEOS-5 GCM 

21 The Fortuna 2.5 version of the GEOS-5 GCM is documented by Molod et al. (2012); it 

22 describes the model performance with several new updates to the earlier MERRA version 

23 (Rinecker et al., 2008). Briefly, it employs Relaxed Arakawa-Schubert scheme (RAS) due to 

24 Moorthi and Suarez (1992) for moist convection with PDFs of cloud water for cloud 

25 microphysics. RAS produces prognostic cloud-cover, and diagnostic ice mixing ratios, and 

26 cloud water. Other upgrades comprise of large-scale condensation and evaporation, auto- 

27 conversion and accretion of cloud water and ice, sedimentation of cloud-ice, and re- 

28 evaporation of falling precipitation following Bacmeister et al. (2006). Its long-wave radiative 

29 transfer calculations are due to Chou and Suarez (1994), and its shortwave radiative transfers 

30 are due to Chou (1998, 1999). These handle interactions with simulated clouds, cloud water, 

3 1 water vapor, and externally prescribed trace gases. In addition, shortwave calculation includes 

32 absorption, scattering and transmission by aerosols, i.e., it treats the direct effect of aerosols 
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1 only. For full details, refer to Molod et al. (2012) and Rinecker et al. (2008). We will refer to 

2 the AGCM as baseline GEOS-5 AGCM. 

3 2.2 McRAS-AC 

4 The latest version of McRAS (Sud and Walker, 2003a) is chosen as the cloud scheme for 

5 including the aerosol-cloud interactions. Its evolution dates back to development and 

6 evaluation of McRAS built from RAS moist convection. McRAS uses cloud microphysics 

7 based on the work of Sundqvist (1988) and Tiedtke (1993) along with other upgrades namely 

8 rain-evaporation (Sud and Molod, 1988) and convective downdrafts (Sud and Walker, 1993). 

9 It is extensively evaluated in a Single Column Model (Sud and Walker, 1999a, Ghan et al., 

10 2000, Xie et al., 2002). Its climate simulations with GEOS-2 GCM (Sud and Walker, 1999b) 

1 1 are more realistic than that of the baseline mode usedl. It produced reasonable intra-seasonal 

12 oscillations (ISO) in GEOS-2 and GEOS-3 AGCMs. The ISOs were well reproduced even in 

13 the NCAR implementation (Maloney and Hartmann, 2000). Nevertheless, these GCM 

14 applications also highlighted some weaknesses that were addressed in subsequent upgrades 

15 (Sud and Walker. 2003a, b). Without aerosol-cloud interaction to provide LPNC or ICNC, 

16 McRAS used empirical equations by Sundqvist (1988) for estimating precipitation production 

17 rate as a function of cloud water, ambient temperature, and cloud type. For radiation, in-cloud 

18 CPNC were assumed for land and ocean following Del Genio et al. (1996), while the volume 

19 and effective radii of CPs were estimated from another set of empirical assumptions (Sud and 

20 Walker, 1999a). 

21 The current aerosol-cloud interaction microphysics modules are documented in Sud and 

22 Lee (2007). The new option is the use of Barahona and Nenes (2009a) as an alternative to Liu 

23 and Penner (2005) aerosol nucleation for IN. The version of McRAS-AC used herein has 

24 McRAS cloud-scheme plus Fountoukis and Nenes (2005) aerosol activation parameterization 

25 to yield CCN and Barahona and Nenes (2009a) scheme for ice nucleation to yield IN and Sud 

26 and Lee (2007) precipitation microphysicsfor liquid and Sundqvist (1988) for precipitation 

27 emerging from mixed phase and ice clouds. Together, these are the sources of LPNC and 

28 IPNC for water and ice clouds, respectively. In-cloud evaporation and/or precipitation and 

29 self collection of cloud water are parameterized by Sud and Lee (2007), which is a recast of 

30 Seifert and Beheng (2001) to obtain relations for thicker clouds encountered in a coarse 

31 resolution GCM. Any change in cloud mass by condensation/deposition and subsequent 

32 removal by precipitation, works interactively through an implicit backward numerical 
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1 solution that is an approximation for the nonlinear coupled differential equations that are 

2 impossible to solve without iteration. Lacking an IPCN based snow precipitation scheme, 

3 McRAS-AC currently uses Sundqvist (1988) parameterization for mixed phase and ice phase 

4 precipitation. However, inclusion of ice nucleation (IN) (Barahona and Nenes, 2009a, b) and 

5 Bergeron-Findeisen cloud water-to-ice mass transfer (Rotstayn, 2000) cloud liquid and ice 

6 mass fractions and corresponding CPNC (LPNC plus IPNC) are calculated consistently to 

7 conserve mass and IPNC budgets. Nevertheless, CPNC reduction is non-linear and is based 

8 on a curve-fitted relationship between cloud mass and number concentration for an assumed 

9 Gamma distribution of cloud particle sizes. Homogenous freezing of in-cloud water drops 

10 surviving until -38°C is enforced through instantaneous freezing of LPNC into IPNC. The 

11 above processes, carried out by several modules, provide an end-to-end treatment of 

12 prognostic cloud water mass (apportioned between liquid and ice) and CPNC, and 

13 precipitation into liquid and snow hydrometeors. Its present implementation is shown through 

14 a block diagram, Figure 1. It shows how aerosol activation by vertical ascent or its equivalent 

15 cooling creates condensate while Sud and Lee (2007) and Sundqvist (1988) cloud micro 

16 physics create precipitation and reduce the CPNC. McRAS-AC is implemented as an option 

17 to work with the baseline GEOS-5 GCM. For more details, the reader mal like to refer to the 

1 8 original papers referenced above. 

19 

20 3. Simulation Experiments 

21 We performed two 10-year long simulations, one with the baseline Fortuna 2.5 version of the 

22 GEOS-5 AGCM (referred to as baseline) with its own cloud scheme and one with the same 

23 GEOS-5 AGCM, but for the newly implemented McRAS-AC cloud physics module as an 

24 option in the AGCM (See details in Table 1 for differences). There are two more differences: 

25 one in the use of simulated clouds without any scaling of cloud fraction, and one is use of 

26 GEOS-5 AGCM PBL cloud treatment instead of the dry convection based in McRAS-AC. 

27 The dry eddy transports pulls up the near surface water vapor to the level of neutral buoyancy 

28 near the PBL top, whereby it raises the height of the cloud base and dries the boundary layer. 

29 In the present implementation, McRAS-AC relies on the PBL scheme of the baseline model. 

30 The monthly climatologies of aerosols are taken from GOCART (Chin et al, 2002) and are 

3 1 based on extensive aerosol model development and calibration/validation exercises at GSFC 

32 (Colarco et al., 2010). Currently we use five externally mixed aerosols namely sulfates, sea- 
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1 salt, mineral dust, black carbon and organic carbon. GOCART model generates time-series of 

2 mass balances of different aerosol species. The aerosol numbers are calculated by making 

3 sectionals of aerosol sizes, called modes. The present goal is to detennine if McRAS-AC can 

4 perform reasonably well in the GEOS-5 AGCM and qualify to perfonn studies of the 

5 influence of aerosols on cloud microphysics and cloud radiative effects (CRE). 

6 

7 4. Results 

8 A comparison of two 10 year integrations, one with McRAS- AC in the GEOS-5 AGCM 

9 (hereafter referred to as “MAC”) and one with the baseline GEOS-5 AGCM (hereafter, 

10 referred to as “CTL”) examine precipitation, clouds, their water paths, and effective radii, and 

1 1 CREs. There are two aspects of this intercomparison; one is the differences due to cloud 

12 parameterization schemes in CTL and MAC, the other is the influence of aerosol activation 

13 and associated cloud particle numbers and sizes; however, the effective radii of cloud liquid 

14 and ice particles are emperically prescribed in CTL runs as functions of temperature and 

1 5 cloud water path. The two most intractable component of the simulated climate change are the 

16 CREs and how these would change with any climate change scenario. GCMs need to simulate 

17 realistic (implying bias free) CREs. The goal here is to detennine how MAC and CTL 

18 climatologies compare with each other and how well they hold up against observations. Can 

19 they simulate the annual cycle reasonably well, and, how much can we trust the model to 

20 produce realistic CRE changes in simulating a climate change scenario of the future? 

21 Specifically the aim is to evaluate MAC seasonal climatology biases and thereafter design 

22 upgrades to ameliorate them. Second, are there reasonable sensitivities to aerosol mass and 

23 number concentration of the real environment, and can they be used to improve model 

24 model’s CREs and understand the influence of aerosols on clouds, and cloud-radiation 

25 interactions and its consequences on the regional climate change? We produced MAC and 

26 CTL seasonal average fields for DJF, MAM, JJA, SON and an annual mean for several key 

27 fields. However, to keep figures-volume low, we only show the climatologies of the two 

28 extreme seasons, DJF and JJA, and the annual means. Although revising the algorithms and 

29 making aerosol input modifications to ameliorate some of the biases are left for Part 2 of the 

30 paper; we will include two test runs that have potential to eliminate both biases of McRAS- 

3 1 AC clouds identified herein. The following are the key highlights of our findings. 

32 


8 



1 4.1 Precipitation Fields 

2 The left two panels of Fig. 2 depict the broad feature of rainfall climatology of MAC and 

3 CTL. Both model-simulations have reasonable ITCZ and SPCZ, with intense convective 

4 rainfall, as expected. In DJF intense rainfall occurs over the South American landmass, 

5 Australia, the tropical islands of Southeast Asia follow the North Pacific Ocean currents, and 

6 the eastern seaboard of continental United States and the Gulf Stream. Large rainfall also 

7 occurs over the rising branch of Ferrell cell between 40S-60S. In JJA, we see tropical rain 

8 including the ITCZ at its northward location. The simulated tropical Pacific ITCZ is 

9 somewhat weaker with less than observed rainfall intensities in the mid-span of the Pacific 

10 ITCZ and somewhat stronger than observed near the land masses at both ends on land, 

1 1 presumably due to orographic intensification of precipitation; thus more water vapor 

12 converges on to land away from its natural location(s) over the ocean where the simulated 

13 rainfall climatology has a deficit. Indian and Asian regions have realistic monsoons and 

14 associated rainfall. Northwards of Sahel, the Sahara desert is dry in JJA as it should be. 

15 Generally, MAC and CTL biases in precipitation are quite similar to each other and hence the 

16 differences from observations are also similar. MAC does somewhat better on the RMSE 

17 scores in DJF and JJA but not on ANN (see Table 3). The majority of the biases are 

18 associated with orographic intensification of precipitation and its related moisture 

19 convergences, a proverbial problem with a number of numerical models that has been largely 

20 solved according to Chao (2012), but the Chao-code is not implemented in the Fortuna 2.5 

21 version of GEOS-5 AGCM. In DJF one sees excessive rainfall biases along southern Andes, 

22 hilly regions of South Africa, and tropical islands of south East Asia. Precipitation biases 

23 around eastern regions of Himalayas, and Andes through South and Central America (red 

24 color in the difference maps). MAC or CTL simulation minus observations, the rainfall 

25 biases are more positive over the tropical Pacific ITCZ in CTL versus MAC whereas the 

26 biases are quite similar in the JJA over the tropical Pacific ITCZ except that lighter colors in 

27 MAC minus OBS means that the correponding biases are lesser in MAC. Overall, both 

28 simulations, MAC and CTL, do reasonably well in comparison to GPCP precipitation data. In 

29 the boreal winter (DJF) season, MAC (CTL) simulate global mean precipitation of 2.89 (2.84) 

30 mm/day versus the somewhat smaller value of 2.68 mm/day in the GPCP data. The 

31 corresponding boreal summer (JJA) values are 3.01 (3.04) mm/day versus 2.71 mm/day in the 

32 GPCP data. Indeed, the simulated precipitation values are consistent with global mean surface 

33 evaporation. Accordingly, global condensation heating of McRAS-AC are 6.0 (8.7) Wm' in 
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1 DJF (JJA) larger. Since SSTs are prescribed, excessive evaporation over the oceans can occur 

2 without any negative feedback that could reduce the SSTs evaporation. (See Table 3 for the 

3 details of data). 

4 DJF averages in the tropics show that the MAC (CTL) rainfall distribution over the sharp 

5 ITCZ is less (more) intense than the GPCP data. However, MAC simulations make up for the 

6 reduction with small increases over the grid cells north and south of the ITCZ. The orographic 

7 rainfall intensification biases are consistently positive and quite similar in both MAC and 

8 CTL simulations. Clearly, the GCM has a problem reproducing observed precipitation with 

9 flow across the steep hills; this is a persistent bias whose solution has remained elusive. For 

10 JJA, both MAC and CTL simulations the Equatorial Pacific ITCZ in the middle of its east- 

1 1 west span is weak similarly, even though CTL simulation has a slightly better organized 

12 ITCZ. On the whole, both cloud schemes show significant biases that need some attention. 

13 In NH winter, precipitable water and precipitation yields are relatively small; accordingly, the 

14 magnitudes of DJF biases are also smaller compared to the wanner southern latitudes. 

15 Conespondingly, precipitation biases over the upslope regions of high orography are also 

16 smaller. McRAS-AC with all the extra degrees of freedom turned out to as realistic as the 

17 baseline GCM. In JJA, the CTL Indian monsoon rainfall is slightly better than MAC, but the 

18 structure of the biases is very similar, while the differences in the biases are relatively small. 

19 Thus comparing precipitation biases in MAC and CTL simulations, deficiencies of McRAS- 

20 AC in some regions despite its better CRE does not warrant keeping McRAS in the 

2 1 development phase of climate simulation studies in which CREa and precipitation changes are 

22 vital for understanding climate change. 

23 Moreover, differences that are statistically significant are important. When this is done with 

24 each model simulation minus observations, the orographically enhanced precipitation biases 

25 overwhelmed the outcome, but when the analysis is done on rainfall differences between the 

26 two models runs, we isolate areas that show the influence of moist physics only (Figure 3). 

27 Large differences are notable over the tropical ITCZ that includes East Pacific in DJA and 

28 mid Pacific in JJA; some differences are significant over the tropical Atlantic and Indian 

29 Ocean in both seasons. In these convergence zones MAC precipitation is less than that of 

30 CTL. Naturally, areas of statistical significance on annual mean fields are even smaller in 

3 1 conformity with our understanding that biases often reduce by averaging not only the annual 

32 cycle but also across space and multi-model ensembles even over shorter time periods. 
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1 Overall, the rainfall differences between MAC and CTL in a 2-tailed student t-test at 95% 

2 significance are small and without much structure except for ice-covered Polar Regions. Thus 

3 both schemes produce similar rainfall fields with very similar biases vis-a-vis GPCP data 

4 barring the biases introduced by orographic intensification which really gives large biases 

5 with respect to GPCP data. 

6 

7 4.2 Cloud Fractions and Water Path 

8 Cloud fractions and in-cloud water particle number concentration are two fields that govern 

9 the most important field, the CREs (Section 4.5.1); hence biases in cloud radiation fields of 

10 the shortwave and longwave radiation can be li nk ed to the biases in cloud fraction and cloud 

1 1 water path and effective radius. The vertical alignment of clouds (Lazaros et ah, 2012) or the 

12 prescribed cloud water-path also influences the CREs. These are discussed at length in the 

13 subsections of Section 4.5 

14 4.2.1 Cloud fractions 

15 We show yearly average zonal mean cloud fractions for the atmospheric column as a whole as 

16 total, that is divided into three: high, middle, and low level values (Figure 4) and their biases 

17 and RMSE with respect to ISCCP data (Rossow and Schiffer, 1999) for DJF, JJA and annual 

18 means (Table 3.1). We examined the simulated total cloud fraction that can be directly 

19 compared with satellite observations. The simulated total cloud fractions of MAC are larger 

20 than those of CTL but lesser than the observed. It zonal distribution tracks the observed quite 

21 well. On this measure, MAC simulations produce lesser biases than the CTL. In comparison, 

22 column averaged cloud fraction biases in CTL are negative globally. Whereas MAC biases 

23 are also negative (except for polar regions), but they are half as much as CTL-biases except 

24 for the tropics. If the simulated cloud fraction is far off, the only way to get reasonable CRE is 

25 to tune cloud water or scaling the cloud fraction and/or their optical properties; while MAC 

26 uses cloud water path scaling, CTL perfonns convective cloud fraction scaling (see Table 1). 

27 Whereas CTL total clouds are consistently less than observed, MAC total clouds are similar 

28 to CTL in the tropics but larger than CTL in the higher latitudes (Figure 4a) and get even 

29 more than the observed in Polar Regions, a bias that has been highlighted before. Further 

30 details of high middle and low level clouds are as follows. 
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1 High clouds are equally robust and realistic in both CTL and MAC simulations except for the 

2 Polar Regions, where MAC high clouds show larger bias than the CTL (Figure 4b). On 

3 seasonal time-scales the biases do look quite different. Some biases of lesser than observed 

4 high cloud over north of India and western United States, and Argentina are common to both 

5 MAC and CTL simulations in DJF as well as JJA (not shown) leading to discernible 

6 signatures in the annual averages. Over the oceans, MAC does better than CTL, which 

7 overestimates high clouds. MAC mid-level clouds are much more realistic than those of CTL 

8 clouds at most places (Figure 4c). The mean RMSE errors of MAC (CTL) for DJF, JJA, ANN 

9 averages are smaller (larger) vis-a-vis ISSCP data as observations (see Table 3). Generally, 

10 CTL produces lesser than observed mid-level clouds. On low level clouds, MAC does a better 

11 job than CTL between 30S to 30N, but at higher latitudes MAC shows huge biases 

12 particularly in the Polar Regions where observations are not so reliable. Regardless, this 

13 makes low level cloud fractions of MAC somewhat inferior. Specifically, large biases in low- 

14 level cloud fractions at high latitudes are seen in all the seasons as well as in the annual 

15 means. CTL simulations produce lesser than observed low-level clouds, but they get better at 

16 higher latitudes (Figure 4d). In summary, clouds above 400hPa, called high clouds, are more 

17 than observed in both MAC and CTL simulations in both DJF and JJA seasons. In DJF, MAC 

18 (CTL) simulates about -0.3 % (0.7%) bias in the observed high clouds fraction of about 

19 22%. In JJA too both MAC (CTL) simulate 0.5 %( 1.3%) more high clouds and those biases 

20 are reasonable in GCM applications. Perhaps RAS, the basic convective scheme of both the 

2 1 CTL and McRAS-AC cloud scheme is responsible for such close similarity in cloud fractions 

22 in the tropics. A word of caution about ISSCP data is in order here. Its retrieval algorithm is 

23 unable to detect very thin cirrus; therefore its bias is towards lesser deep clouds and could 

24 explain high cloud biases as a spurious model deficiency. 

25 The differences in the middle level, 700-400hPa, clouds are very much better in MAC vis-a- 

26 vis CTL. MAC biases are 2.2% for DJF and 1.5% in JJA versus -10.3% for DJF and -10.3% 

27 for JJA for CTL simulations. MAC scheme started producing more mid-level clouds after 

28 melting of snowfall at 0°C was introduced (Sud and Walker, 2003b). The snow melts around 

29 500hPa in the tropics, produces an inversion that debars cumulus towers from penetrating 

30 through it. We also like to point out that mid-level cloud fraction in the ISSCP data may be 

31 too large (Chen and DelGenio, 2008) because the mid-level clouds percentages are almost 

32 same as high clouds percentages, but for now MAC simulates them and our radiation balances 
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1 suggest that cumulus congestus and mid level detrainment by high latitude cumulus clouds 

2 (Johnson et ah, 1999) may be reasonable; however, the fall velocity of an embrionic 

3 hydrometeors is chosen for optimal estimates of cloud water (Sud and Lee, 2007). Even 

4 though both large-scale schemes use Slingo and Ritter (1987) type of critical relative 

5 humidity dependence for onset of stratiform clouds, the McRAS-AC determines only its 

6 tendency and not the amount (Sud and Walker, 1999a) as well as employs moist convection 

7 associated with stratiform clouds as an additional upgrade (Sud and Walker, 2003b). The 

8 latter would tend to increase the cloud fraction and its vertical correlation that shows up in the 

9 CRE analysis of Oreopoulos et al. (2012). MAC simulates almost 50% more low level clouds 

10 and those are all in high latitude regions. Similar high latitude biases in the CTL case are half 

1 1 as much. Overall, CTL simulated PBL clouds are about 4% less than the observed. However, 

12 in the tropics, MAC simulation is still as good. McRAS-AC module in all the past 

13 applications had a PBL cloud scheme that was based on dry convection and was linked to 

14 surface fluxes. In GEOS-5 GCM application used in MAC simulation, we rely on the PBL 

15 cloud scheme of the baseline GEOS-5 GCM. There seems to be need for reworking on the 

16 PBL clouds. In the regions with moist convection, moisture exits near the top of the PBL 

17 through dry convection but at high latitude with little or no convection, the moisture saturates 

18 the PBL and makes excessive PBL clouds. For example, in the Polar Regions, often 

19 precipitation emerges as tiny snow particles from the ice fog falling out of clear sky (diamond 

20 dust; Greenler, 1999). Our cloud model, McRAS-AC, may be identifying it as fully cloudy 

2 1 because local RH of the ambient atmosphere exceeds the saturation vapor pressure for ice, the 

22 criteria used to identify a large scale cloud in the model. Moreover PBL clouds develop over a 

23 long period of time and that gives them more time to entrain the aerosols into the turbulent 

24 plume. These are some of the ideas that can mitigate the high latitude cloud biases of McRAS 

25 that reflect in MAC simulations. 

26 

27 4.2.2 Total Cloud Water Path 

28 The geographical distribution of total cloud water path (= sum of liquid plus ice water paths) 

29 simulated by MAC and CTL are shown in Figure 5. First, CTL biases are consistently 

30 negative in the high latitude regions and positive in the tropics. Most of the positive moisture 

3 1 biases are in the tropics, but still it is not able to counter the high latitude negative biases (see 

32 Table 3 for means and RMSEs). However, since the effective radii of cloud drops are 
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1 prescribed in CTL, it uncouples the cloud mass and number density from the cloud water 

2 budget and allows the modeller to optimize the results bu tuning. Indeed the MAC simulation 

3 too has its biases. Its RMSEs are similar to those of CTL whereas its biases are better 

4 connected with the cloud water path. Too high total and liquid cloud water path is due to high 

5 water content in the storm tracks along the North Pacific Currents and the Gulf Stream over 

6 the Eastern boundaries of the Asian and North American continents. This bias is also reflected 

7 in the annual averages. The DJF biases are much smaller although the RMSE are not. This is 

8 mostly related to negative biases in Polar Regions, less than observed cloud water path in 

9 convective regions, lack of stratus clouds of the west coast of North and South America. 

10 

1 1 4.2.3 Liquid Cloud Water Path 

12 The simulated liquid cloud water path is shown in Figure 6. Here OBS stands for SSM/I 

13 dataset from Greenwald et al.( 1993) and Weng and Groudy,1996 which is good for oceans 

14 only. This has a distinctly different pattern than the total water path. Its mean statistics are 

15 better for MAC than CTL, while its RMSEs are not. Evidently, MAC produces less cloud 

16 water in the ITCZ, SPCZ in DJF, JJA and ANN averages; in other words, McRAS-AC 

17 generates more than observed precipitation in convective areas. Once that happens, both mass 

18 and number densities of cloud water reduce proportionally. Moreover, our analysis of 

19 fractional cloudiness shows that low level cloud fractions are more than the observed; 

20 therefore large water content in the stratifonn regime and smaller particle sizes may also have 

21 contributed to this bias. However, adjusting the tuneable parameters of the scheme can 

22 potentially ameliorate this problem. On other specifics, in DJF there is too much liquid as 

23 well as total water in the 40-60S region of roaring winds. In fact, as the total cloud water bias 

24 reverses in JJA while the liquid cloud water still has a positive bias, we infer that the ice 

25 amount over the region in JJA is too small and that points to lack of IN and a delayed 

26 Bergeron process which waits to kick in until sufficient IN are available for water vapour to 

27 deposit on. This may well be related to lack of ice nucleating aerosols. Since ice is also an 

28 absorber of solar radiation, lack of ice in the clouds near the boundary layer may be one of the 

29 causes for the bias in solar absorption, which can only appear in DJF because the solar 

30 radiation at 40-60S latitudes is much smaller in JJA. Low cloud water path off the west coast 

31 of Americas in JJA is expected lacking the boundary layer stratus and that indeed happens in 

32 both MAC and CTL simulations. The liquid water content in the storm track regions from 
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1 Asia to west coast of North America (over the ocean) has a strong positive (negative) bias for 

2 MAC (CTL). In fact if we examine the total water bias over the Gulf Stream and North 

3 Pacific Current, the large positive bias in liquid and total cloud water in MAC simulation is 

4 notable. This is related to large-scale clouds and presumably high aerosol content of the 

5 ambient air mass that make precipitation rate too low. 

6 

7 4.2.4 Sensitivity of Cloud Water to Aerosols 

8 A one year test simulation is made with interactive aerosols generate by the GOCART model 

9 instead of the climatologically prescribed aerosols generated by the same model. The 

10 interactive run reduces the aerosols within the storm in the stonn track regions. Rain storms 

1 1 not only remove the activated aerosols as CCN and/or IN through precipitation, but also 

12 scavenge the inactivated ones through wash down, which is difficult to parameterize. In the 

13 one year test, the high cloudiness and reflectivity bias of the clouds in the storm track regions 

14 are greatly reduced giving us hope that the problem can be solved without a major research 

15 effort (Figure 7). 

16 

17 4.3 Cloud Particle Effective Radii 

18 The effective radii of ice and liquid cloud particles of MAC and CTL simulations are 

19 compared against satellite data (Platnick et ah, 2003). It must be pointed out that the 

20 validation satellite data based on radiances is unable to see the entire vertical distribution of 

21 cloud particles and their effective radii and therefore these radii are based on what the 

22 instrument sees near the top. Figure 8 shows MAC and CTL effective radii for liquid and ice 

23 cloud particles. The liquid effective radii of MAC simulation are in the range of 10-18 pm 

24 with a global average of 14 to 14.4 pm in DJF and JJA seasons. It is too small by 1-1.5 pm, 

25 which is well within the spread of the observations; inferring the same from empirical 

26 relations used in CTL, the values are in the 9-14 p m range with global average value of 10. 1 

27 to 10.5 pm for DJF and JJA respectively. Thus the CTL effective radii are about 30% smaller. 

28 Indeed MAC is also biased somewhat on the smaller size, but the size depends upon the 

29 aerosols activation and we believe, the input aerosols easily may have more than 5-10% 

30 biases globally because their numbers are not well validated over the oceans. The fact that we 

3 1 get liquid effective radius so close to observation shows that on this measure liquid aerosol 
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1 parameterization of MAC is as good as it can be expected of a realistic cloud model. Because 

2 of the assumptions about the sub-grid scales, and having seen the biases in clouds, the 

3 simulated zonal averages effective radii are as reasonable as one can hope to get in a GCM, 

4 whereas the prescribed effective radii in CTL runs are biased much more. 

5 The simulated effective radii of cloud ice particles 24-42 pm with an average 28.3 to 29.9 pm 

6 for JJA and DJF respectively. Clearly MAC simulates larger effective radii but again it is 

7 about 10% or 3 pm too large and that too is well within the observational spread. In fact 3 pm 

8 in 30 pm can be easily explained by the range of variations between the effective and volume 

9 radius of ice particles. Again in CTL, the empirically estimated ice particle effective radius is 

10 ~3 p m less than the observed, however, its zonal biases tell us that polar ice clouds have 

1 1 problem and that is consistent with their low number and large particle size. The RMSE errors 

12 of MAC are twice as large but that is to be expected because MAC numbers, being actual 

13 predictions, have a large spread whereas CTL values do not; (see Figure 8). Based on these 

14 results, MAC is simulating liquid and ice cloud effective radii reasonably well and is ready to 

15 perform as an aerosol-cloud-radiation interaction model with its current cloud physics. The 

16 grid average liquid and ice particle numbers are about 40 and 3.8 cm' with corresponding in- 

17 cloud values of (—90 and 10 cm' ). There are no observed values for these but judging by the 

18 radiation imbalances at the top of the atmosphere, we assume that McRAS-AC simulated 

19 numbers have large zonal biases that may be aerosol related. 

20 

21 4.5.1 TOA radiation budget 

22 In this subsection we assess briefly the verisimilitude of the radiation budget produced by the 

23 10-yr simulations of GEOS-5 with the two cloud schemes. We compare model simulated 

24 TOA LW and SW zonal fluxes to their counterparts from CERES products (Loeb et al. 2009, 

25 CERES data set EBAF 2.6). Figure 9 shows for DJF, JJA, and ANN fields while global mean 

26 biases and RMS errors are displayed in Table 3. Assessments of deficiencies in the simulated 

27 cloud fields that lead to discrepancies between the simulated and the observed radiative fields 

28 are mostly left for the next subsection which frames the discussion in terms of cloud radiative 

29 effects (CRE). 

30 The zonal-average radiation plots shown in Figure 9 indicate that for the most part, MAC 

3 1 matches OBS better than CTL. As will be discussed further later, while its tropical convective 

32 clouds have too little cloud water/ice and somewhat lower height, the resulting outgoing 
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1 longwave radiation (OLR) and absorbed shortwave radiation (ASW) are closer to 

2 observations than the CTL whose overactive convection yields too much reflected SW and 

3 too little OLR. MAC is less skillful in midlatitudes, however, with too much ASW in the SH 

4 summer (discussed extensively below) and too much OLR in the NH summer. 

5 Overall, both schemes are in reasonable agreement with observations in terms of global mean 

6 OLR (Table 3), but MAC RMSE values are better by ~2 Wm' for both seasonal and annual 

7 averages, indicating larger spatial error cancellations for CTL. The annual averages of global 

8 ASW are about the same for both cloud schemes, and about 2.5 -3.0 Wm' larger than the 

9 observations. However, summer and winter global values differ substantially between the two 

10 schemes. MAC’S ASW differs by ~1 7 Wm" between DJF and JJA while the corresponding 

1 1 differences in both observations and CTL are about half. Again, this is a result of deficient 

12 simulation of SH midlatitude marine clouds by MAC (the reasons for this bias and possible 

13 solutions are discussed later) making the RMSE slightly worse than CTL in DJF (Table 3). 

14 The same deficiency impacts also the global net TOA flux which is — 9 Wm' too high for 

15 MAC in DJF compared to the CERES. The global net TOA flux for JJA is within ~1 Wm' 2 of 

16 observations for both MAC and CTL, but the DJF error of MAC is too large, yielding a 

17 substantial excess of 8.5 Wm' in annual global net TOA radiation. Our simulations are not 

18 much affected by this large energy imbalance because of externally prescribed SSTs. The 

19 RMSEs of net TOA radiation are worse for MAC than CTL for DJF, but are better for JJA, 

20 and about the same for ANN . Nevertheless, taking all radiation quantities into account, and 

21 focusing on their RMSEs (which are not affected by spatial cancellations), MAC generally 

22 produces radiation fields that are closer to observations than CTL. This is reaffirmed in the 

23 next subsection which examines cloud radiative effects. 

24 

25 4.5.2 TOA Cloud Radiative Effect 

26 A well-established way for assessing the influence of clouds on the radiation budget is via the 

27 cloud radiative effect (CRE), a quantity also popularly known as cloud radiative forcing 

28 (Harrison et al. 1989). CRE for either solar/shortwave (SW) or thennal infrared/longwave 

29 (LW) radiation is defined as: 

Tjcld J7 clr 

30 LW,SW r LW,SW r LW,SW (la) 

3 1 which can also be recast as 
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1 CRE lwsw = C lol (F- sw -F^ sw ) (lb) 

2 under the assumption that the cloudy sky flux can be written as the linear combination of clear 

3 and overcast fluxes. In the above, F is the net downward (i.e., downward minus upward) flux 

4 (LW or SW), the superscripts clr designates clear (cloud-free) skies, eld designates all-sky 

5 conditions (containing a mixture of cloudy and clear skies), and ovc designates overcast skies 

6 (100% cloud fraction); C tot is the total vertically projected cloud fraction which in the AGCM 

7 depends on individual layer cloud fractions and the cloud overlap assumptions. While both 

8 definitions can be used for analyzing observational data, the model CRE always comes from 

9 eq. (la). Nevertheless, eq. (lb) is preferable for interpreting AGCM CRE. For two different 

10 cloud schemes producing the same C tot , the CRE differences mainly arise from their water 

1 1 path/effective radius differences (their combined effect is captured by the cloud optical depth) 

12 in the SW, and cloud top height differences in the LW (although optical depth differences also 

13 play some role at low values of optical depth) through their effect on F° vc . 

14 The CRE as defined above can be calculated at either TOA or at the surface. Here we only 

15 show TOA results for which more reliable observed values are available. In the SW, the CRE 

16 TOA is usually negative because the net (absorbed) flux for cloudy skies is smaller than for 

17 clear skies. In the LW, the TOA CRE is usually positive because the upward TOA flux is 

18 greater under clear skies than cloudy skies (the downward flux is zero in both cases). The 

19 CRE of net radiation is, CRE net = CRE sw + CRE LW , and can be positive or negative 

20 depending on cloud type. Measurements of TOA CRE are readily available from CERES, 

21 among other sources, and can be used for model evaluation. We use the EBAF v. 2.6 of the 

22 CERES data set (Loeb et al. 2009). 

23 

24 In an earlier paper, Oreopoulos et al. (2012) showed that for diagnostic radiation calculations 

25 with a different radiation scheme, the TOA CRE and its sensitivity to the vertical distribution 

26 of clouds was very different for the CTL and MAC clouds. Based on results of Oreopoulos et 

27 al. (2012) we naturally expected substantial CRE differences between MAC and CTL. The 

28 TOA CRElw differences between observations and the two model runs are shown in Figure 

29 10. In this plot as well as similar ones that follow for CREsw and CRE net , red and blue colors 

30 represent positive and negative biases with color intensity proportional to the magnitude of 

31 the bias. The lighter shade of colors for MAC minus OBS compared to CTL minus OBS is 
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1 indicative of smaller biases for MAC, and this is especially true over the Pacific. The CRE LW 

2 biases in both simulations are consistent with those of cloud fraction bias discussed earlier in 

3 subsection 4.2.1. In the ITCZ region, MAC exhibits smaller CRElw than observations, 

4 probably because its convective clouds are too low or too thin, while CTL exhibits the 

5 opposite behavior, i.e., larger than observed CRE LW , suggesting that convective clouds in 

6 CTL may be too high and/or too thick and too spatially extensive. In general, the MAC 

7 underestimates are lower than the CTL overestimates. Lower cloud tops in MAC may be due 

8 to the influence of quadratic entrainment in McRAS-AC (and neglect of convective heigh 

9 increase by freezing of cloud ice and precipitating hydrometeors) versus linear in the standard 

10 RAS of the GEOS-5 GCM (Sud and Walker, 2003a). Larger entrainment aloft will reduce the 

1 1 in-cloud moist static energy and will keep them shallow. In fact, when we examine all the 

12 places with abundance of convective clouds, the simulated CRE LW in CTL is consistently too 

13 large suggesting that cloud tops are too high or the cloud free areas are too few or too small. 

14 Because lesser entrainment aloft will need lesser cloud mass to neutralize the cloud work 

15 function above the critical value, which would reduce the cloud fraction and mitigate the 

16 lesser CRE LW biases as opposed to making them worse. Despite the previously mentioned 

17 weaknesses in the simulation of southern midlatitude ocean clouds by MAC, CRElw biases 

18 are not as high because most of the clouds at these latitudes reside in the lower troposphere 

19 and do not have much influence on the CRElw. On the other hand, biases in the snow and ice 

20 covered polar regions (where observed CREs may be less reliable), both the positive and 

21 negative biases are generally larger in MAC than in CTL for reasons that remain currently 

22 unidentified. 

23 Overall, the CTL scheme underestimates global CRElw by ~4 Wm' (Table 3) mainly because 

24 of the systematic underestimates in convective areas. The MAC simulation approaches the 

25 global observed value within ~ 1.0- 1.5 Wm' and achieves better RMSE scores than CTL in 

26 both the seasonal and annual means by — 1.5-3 Wm' . Figure 1 1 shows TOA CREsw radiation 

27 difference maps. CREsw fields of CTL minus OBS have deeper colors with more structure 

28 compared to those of MAC minus OBS. Large differences in the biases are evident in MAC 

29 and CTL over northern midlatitudes in JJA. But the most prominent MAC biases 

30 (underestimates) appear in DJF within the 40S-60S latitude zone where MAC produces too 

31 few CPNC that have consequently too large sizes (see Section 4.3.). Since the region is 

32 dominated by sea salt aerosols, we hypothesize that either these aerosols are not activated 

33 adequately, or the inferred particle numbers from the GOCART mass concentrations are too 
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1 low. To examine the impact of the latter possibility, we conducted one year run where we 

2 reduced the sea-salt aerosol diameter by 50% across the board resulting in an 8-fold increase 

3 in aerosol particle number density (APND). This is a reasonable test because GOCART 

4 simulates mass balances employing only the mass tendency as sum of sources, sinks, aerosol 

5 chemistry and advection; APND for different bins is estimated from volume radius and 

6 density that match the aerosol optical thickness. The 8-fold increase in the sea salt APND 

7 resulted in a TOA CREsw field very similar to that of the CTL (Figure 12). While this 

8 experiment isolated the cause of the bias, it cannot be considered a solution. Greater ice 

9 particle numbers can also be created by a physically based ice-cloud particle splintering 

10 algorithms. The region is predominantly in the rising branch of the Ferrell cell where winds 

11 are strong and gusty, consequently CPNC increases due to cloud particle colliding and 

12 shattering, ignored in the current version of McRAS-AC, can be significant. Another 

13 mechanism that would increase IPNC is liquid cloud particles glaciating sooner as opposed to 

14 being depleted by Bergeron-Findeison mass exchange between water drops and ice particles 

15 through evaporation-desublimation process. Eliminating the biases with better algorithms, 

16 would not only mitigate the CPNC biases over 40S-60S, but would have the potential benefit 

17 of eliminating CREsw biases elsewhere as well. We are actively working on a physically- 

1 8 based solution to this problem. 

19 The similarity of some biases appearing in both simulations suggest either the influence of 

20 their common RAS (Moorthi and Suarez, 1981) heritance or other shared model deficiencies 

21 such as absence of boundary-layer stratus clouds, and excessive orographic precipitation. 

22 Wherever the diagnostics show significantly similar biases in MAC minus OBS and CTL 

23 minus OBS, a common cause, not related to aerosol-cloud interaction, is assumed to be the 

24 culprit. Regarding the positive biases (underestimates of CRE S w) over the PBL stratus regions 

25 off the west coast of California and Peru, an ad hoc tuning of the PBL moisture transport in 

26 the vertical ameliorated this problem in a test version of GEOS-5 GCM, but a more physically 

27 sound alternative is needed. Such is the scheme of Bretherton and Park (2008) which was 

28 successful in simulating realistic stratus clouds off west coast of north and south America as 

29 shown in Kay et al. (2012). The current baseline GEOS-5 GCM lacks PBL stratus, and in this 

30 exercise, both MAC and CTL simulations exhibit similar CREsw biases due to this inherent 

3 1 flaw(s) in the model’s PBL convection. 

32 
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1 The positive CREsw biases in southern midlatitude oceans beneath the Ferrell cell between 

2 40S-60S are big enough to cause a global mean underestimate of 5 Wm' in DJF when the SH 

3 insolation peaks. For the same reason, for SH summer, MAC’S RMSE is slightly larger than 

4 CTL even though for JJA and ANN the RMSE is notably smaller for MAC. The global ANN 

5 CREsw of MAC is 1 Wm' too low since the SH summer underestimate is larger than the NH 

6 summer overestimate. CTL simulates better the summer SH, but has in general more bias 

7 compensations as evidenced by the larger RMSEs in JJA and ANN. 

8 

9 The bias fields of CRE net (Figure 13) reflect previously discussed issues: In areas where 

10 CRElw is small, the CREsw biases take over, see for example the SH midlatitude oceans 

1 1 (MAC) and PBL stratus areas (both simulations). MAC fares better in the intensely 

12 convective regions: apparently its CRE S w and CRE LW underestimates largely cancel out 

13 because they have opposite signs. On the other hand for the CTL the tropical overestimates of 

14 CREsw are significantly larger than the overestimates of CRElw resulting in too strong (too 

15 negative) CRE net , thus implying that the region loses radiative energy at a rate larger than that 

16 of CERES observations. Based on the global values of CRE net alone (Table 3), one would be 

17 tempted to conclude that CTL simulates better cloud fields than MAC. But much of the 

18 agreement with CERES is fortuitous and results from cancellations between the SW and LW 

19 CREs as well as spatial cancellations. Indeed, the MAC RMSEs of CRE net are lower on both 

20 seasonal and annual basis. 

21 

22 4.6 Comments on the Statistics of Circulation 

23 At the outset, one notes that the MAC and CTL simulations are closer to each other than the 

24 Satellite data “so called observations” each of which provides useful guidance for model’s 

25 biases but each data has its own uncertainties. We first like to detennine where McRAS-AC 

26 really made a difference and whether it is statistically significant and/or beneficial for the 

27 climate forecast. We found some differences in the precipitation and circulation, but most of 

28 them were quite local and without much large scale structure. Other differences were in 

29 regions where the input data to the 4DDA analysis system are sparse, the analysis largely 

30 supported the background model simulations; clearly, there the biases were lesser for CTL as 

31 opposed to MAC because analysis used the baseline GEOS-5 GCM. Its example is 100 to 200 

32 hPa temperatures. We decided to postpone this analysis for Part 2 of this paper or until some 
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large biases of McRAS-AC are removed. This will figure in our future research. Significant 
changes in precipitation in convective regions are the only easily interpretable differences 
between MAC and CTL and those were discussed in Section 4. 1 .2 


5 Summary and Conclusion 

We examined 10 year long simulations with the GEOS-5 GCM with prescribed SSTs from 
SST analysis (reference). One used the baseline model and one used the McRAS-AC cloud 
physics, which includes McRAS cloud physics developed by Sud and Walker (1999a) with 
follow-on upgrades (Sud and Walker, 2003a, b) to include aerosol-cloud-radiation interactions 
(Sud and Lee, 2007). The New McRAS-AC has Barahona and Nenes (2009a, b) ice nucleation 
that is substituted for Liu and Penner (2005) ice nucleation scheme. Results show: 

a) McRAS-AC (MAC) simulation produced as good circulation and precipitation fields as 
the baseline GEOS-5 AGCM (CTL) simulation. There are patchy areas of significant 
differences in the circulation and precipitation, but most of the major circulation features 
are very similar. Accordingly, it is difficult to say categorically which one is better. In the 
mean and RMSE biases of precipitation, MAC simulation is a little better. Nevertheless, 
large 40S-60S biases in radiative CREs and cloud water path over the stonn track regions 
and missing low level stratus does not justify declaring MAC to be better than CTL. 
However, for addressing these issues, we have identified good strategies that are 
supported by the subsequent sensitivity tests. Specifically insufficient cloud particle 
numbers over the 40S-60S regions are related to extremely deficient sea salt aerosols as 
well as lacking cloud particle enhancement by collision and splintering. While the high 
cloud water path in the storm track region is related to inefficient wash-down of aerosols 
by precipitation over the storm track regions of the North Pacific Currents and the Gulf 
Stream. 

b) MAC simulated better surface and TOA short and longwave radiative fields and their 
CREs but again its potential benefits were largely mitigated by large biases in the 
aforementioned regions. There are several plausible options for eliminating these biases, 
but we need to detennine, which of them are most defensible physically. Elsewhere 
MAC produced better cloud fractions, cloud and ice water paths. 

c) The original McRAS-AC application produced shallow clouds based on dry convection, 
but in the GEOS-5 GCM application it was dropped in favour of the current GEOS-5 
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1 GCM PBL physics. This compromise may be responsible for McRAS-AC simulating so 

2 poor PBL clouds in high latitudes and keeping the low level clouds close to the surface. 

3 d) The cloud particle effective radii for water (ice) clouds is reasonable (somewhat larger) 

4 as compared to the observations, but it is definitely better than the empirical estimates of 

5 the baseline GEOS-5 GCM. Cloud ice particle splintering has the potential to ameliorate 

6 this bias. 

7 Recognizing that aerosol cloud radiation interaction parameterization is in its infancy and 

8 most of the present day models still have discernible biases in aerosol cloud interactions (Kay 

9 et ah, 2012), we submit that McRAS-AC cloud scheme is able to perform as well as the 

10 baseline cloud scheme of the GEOS-5 GCM despite a few large regional biases that are 

1 1 potentially correctable but reduces its RMSEs of the key fields. At this time, the perfonnance 

12 of McRAS-AC in GEOS-5 GCM is as good as the baseline GEOS-5 GCM, if not 

13 unequivocally better yet, consequently one can start using McRAS-AC modified GEOS-5 

14 GCM for climate research in which the focus is on studying the CRE of clouds as opposed to 

15 waiting for next upgraded that can make McRAS-AC more realistic. 

16 Biased CREs are the primary source of biased background circulation in a climate simulation. 

17 Sometimes the climate change problem has smaller radiative forcing anomalies than the 

18 model’s biases, e.g., doubled CO 2 effects produce a mere 3-4 Wm' radiative forcing 

19 anomaly, while models simulating its influences may have biases as large as 30 Wm' in 

20 many regions. Realizing that even the juxtapositions of model-biases vary among the models, 

21 consensus among their climate predictions is unlikely. Increasing the size of the domain 

22 would help because the CRE biases often reduce in averaging and when the region becomes 

23 sufficiently large, the climate change causing radiative forcing anomalies may exceed the 

24 model's biases and hence the possibility of consensus increases. Its prime example is better 

25 consensus among models in the global means. But if our goal is to capture climate 

26 prediction at the regional scales, we must reduce the CRE biases as much as possible. 
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1 Table 1. Parameterizations in GEOS-5 GCM and McRAS-AC Cloud Scheme (s) 


Model Specifications 

Baseline GEOS-5 GCM 

McRAS-AC in GEOS-5 GCM 

Deep Convection 

RAS (Moorthi and Suarez (1992) 

McRAS (Sud and Walker, 1999a) 

Stratifonn cloud 

Smith (1990), Molod (2012) 

Sud and Walker (1999a) 

Precipitation liquid 

Rienecker et al. (2008) 

Sud and Lee (2007) 

Precipitation snow 

Rienecker et al. (2008) 

Sundqvist (1988) 

Rain evaporation 

Bacmeister et al. (2006) 

Sud and Molod (1988) 

Cloud microphysics 

Single moment microphysics 

Double moment microphysics 

Cloud Scaling 

Convective cloud fraction scaled 
(Rinecker et al. 2008). 

Cloud water Path scaled 
(Sud and Walker, 1999a) 

Aerosol Effects 

Direct effects only 

Both direct and indirect effects. 

CCN Activation 

Not included 

Fountoukis and Nenes (2005) 

IN Nucleation 

Not included 

Barahona and Nenes (2009a) 

Liquid Particle Number 

CP Effectiv radius determined as a 

LPNC tendency budget of sources 

Concentration (LPNC) 

function of T and P 

and Sinks. 

Ice Particle Number 

Not included; Effectiv radius 

IPNC tendency budget of sources 

Concentration (IPNC) 

determined as a function of T and P 

and Sinks. 

Cloud liquid mass 

Emperical equation for cloud water 
nass fraction as a function of T 

Solution of cloud liquid mass 
Tendency Eq. minus BF loss. 

Cloud Ice mass 

Total cloud mass minus liquid mass 

Solution of cloud liquid mass 
Tendency Eq. minus BF loss. 


2 

3 

4 

5 

6 
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1 Table 2. Descriptions of Simulation Experiments 


Simulation 

Descriptions 

Years 

Experiments 



Control Run (CTL) 

GEOS-5 GCM standard moist physics 

10 

McRAS-AC Run (MAC) 

GEOS-5 GCM with McRAS-AC 

10 

MAC run with half the 

McRAS-AC with half size sea salt (eight 

1 year Test 

sea salt particle size 

times number density or 8xSS) 


Interactive GOCART 

GOCART interactive with GEOS-5 

lyear Test 

with MAC 

McRAS-AC 



2 

3 

4 Table 3. Global mean MAC and CTL Simulation Fields versus Observations 


Periods (left to right) 

*OBS:Observed Mean 

MAC: Mean/(RMSE) 

CTL: Mean/(RMSE) 

Fields (top to bottom) 

DJF 

JJA 

ANN 

DJF 

JJA 

ANN 

DJF 

JJA 

ANN 

Precipitation 

2.68 

2.71 

2.68 

2.89 

3.01 

2.92 

2.84 

3.05 

2.89 

(mm day" 1 ) 




(1.54) 

(1.87) 

(1.32) 

(1.58) 

(1.92) 

(1.29) 

Total Cloud Fraction 

67.0 

65.5 

66.4 

56.7 

55.0 

55.6 

44.3 

44.8 

44.5 





(17.1) 

(16.6) 

(15.8) 

(25.6) 

(24.8) 

(24.3) 

High Cloud Fraction 

21.9 

21.9 

21.7 

21.6 

22.4 

22.3 

22.6 

23.2 

23.1 





(8.9) 

(9.3) 

(6.6) 

(8.9) 

(9.8) 

(7.6) 

Middle Cloud Fraction 

19.6 

17.7 

19.2 

21.9 

19.2 

20.5 

9.9 

9.9 

9.8 





(9.3) 

(7.4) 

(6.7) 

(13.3) 

(10.8) 

(11.5) 

Low Cloud Fraction 

24.9 

26.6 

25.6 

37.0 

35.4 

35.9 

22.0 

22.8 

22.2 





(21.2) 

(18.5) 

(19.8) 

(13.7) 

(17.4) 

(13.1) 

Cloud Liquid Water Path 

84.3 

85.8 

84.3 

76.6 

84.2 

79.4 

72.9 

74.5 

72.4 

(g m" 2 ) 




(40.1) 

(44.9) 

(34.1) 

(26.4) 

(33.2) 

(24.2) 

Cloud Total Water Path 

89.9 

90.5 

88.2 

92.1 

107.0 

98.3 

77.2 

82.4 

77.8 

(g m" 2 ) 




(51.4) 

(65.1) 

(41.1) 

(59/8) 

(57.9) 

(48.7) 
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^ ^ l/i 4^ W K) 


Cloud-ice effective Radius 

24.8 

25.6 

25.2 

29.9 

28.3 

28.6 

21.5 

21.9 

21.6 

(pm) 




(8.8) 

(9.3) 

(7.1) 

(4.5) 

(6.8) 

(4.4) 

Cloud-drop effective Radius 

15.2 

16.3 

15.6 

14.3 

14.4 

14.0 

10.1 

10.5 

10.3 

(pm) 




(4.2) 

(4.4) 

(3.3) 

(6.1) 

(7.1) 

(6.1) 

Grid Average/In-cloud 




4.1/ 

3.5/ 

4.1/ 




IPNC (# cm" 3 ) 




10.6 

9.4 

10.7 




Grid Average/In-cloud 




35.0/ 

44.5/ 

44.3/ 




LPNC (# cm" 3 ) 




68.9 

93.1 

90.3 




OLR 

236.9 

243.3 

239.7 

236.0 

242.0 

238.6 

237.4 

245.8 

241.1 

(W m" 2 ) 




(8.7) 

(9.8) 

(7.0) 

(10.3) 

(12.1) 

(9.0) 

ASW 

244.5 

235.7 

240.5 

252.2 

235.3 

243.3 

246.7 

239.0 

243.0 

(W m" 2 ) 




(18.9) 

(15.5) 

(12.1) 

(17.2) 

(21.7) 

(15.6) 

Total TOA Rad. 

7.6 

-7.6 

0.83 

16.2 

-6.6 

4.7 

9.3 

-6.8 

2.0 

(W m" 2 ) 




(17.5) 

(12.4) 

(11.1) 

(13.4) 

(16.3) 

(H.2) 

LW TOA CRE 

25.9 

26.3 

26.2 

24.5 

25.4 

25.3 

21.6 

22.2 

22.2 

(W m" 2 ) 




(7.7) 

(7.4) 

(6.0) 

(9.3) 

(10.2) 

(8.3) 

SW TOA CRE 

-51.6 

-44.8 

-47.3 

-45.6 

-46.7 

-46.3 

-50.8 

-43.2 

-46.4 

(W m" 2 ) 




(17.9) 

(16.2) 

(12.2) 

(17.3) 

(20.9) 

(15.3) 

Total CRE TOA 

-25.6 

-18.4 

-21.1 

-21.1 

-21.3 

-21.1 

-27.8 

-17.2 

-21.7 

(W m" 2 ) 




(16.9) 

(14.2) 

(11.3) 

(18.4) 

(21.0) 

(15.2) 


1 * Datasets deployed 

a) GPCP for Precipitation (Adler et al., 2003) 

b) 1SCCP for Clouds (Rossow, and Schiffer, 1999) 

c) SSM/1 for liquid water path (Greenwald et al., 1993) 

d) MOD1S for effective radii and total water path (Platnick et al., 2003) 

e) CERES for TOA Radiation (Loeb et al., 2009) 

f) SRB for Surface Radiation (Wilber et al., 2006) 
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2 Figure Captions 

3 Figure 1.: Cloud Physics and Microphysics of baseline GEOS-5 GCM and its modification to 

4 allow McRAS-AC to replace the relevant moist physics module of the GCM 

5 Figure 2.: Simulated ten year mean precipitation (mm day' 1 ) for DJF(top) , JJA (middle), and 

6 Annual mean (bottom) in MAC and CTL runs (left 2 panels) and MAC minus OBS and CTL 

7 minus OBS (right two panels) ; GPCP data represents OBS. 

8 Figure 3.: Statistically significant precipitation differences (mm day-1) are color coded: MAC 

9 minus CTL for DJF (top) , JJA (middle), and Annual mean (bottom) using a 2-tailed student t- 

10 test. 

11 Figure 4.: Ten year mean zonal average cloud Fractions for a) entire column atmosphere, b) 

12 400hPa to top, c) 400-700 hPa, d) and 700hPA to surface in MAC and CTL simulations 

13 versus ISSCP and MODIS Climatology. 

14 Figure 5.: Total Cloud Water Path for the full depth of the atmosphere in g m' 2 for DJF (top) , 

15 JJA (middle), and Annual mean (bottom) in MAC and CTL runs and MAC minus OBS (left 

16 two panels) and CTL minus OBS (right two panels); MODIS data represents OBS. 

17 Figure 6.: Same as Figure 5 except for Cloud Liquid Water Path; SSMI data represents OBS 

18 Figure 7.: Total Water Path for JJA in g m' 2 ; interactive aerosol chemistry minus GOCART. 

19 Figure 8.: Zonal Average plots of a) ice and b) liquid cloud particle effective radius (pm). 

20 Different colors are used for Land, Ocean, and Total (colors legends are displayed in the 

21 Plots) for MODIS data (dotted), MAC simulations data (solid), and CTL simulations (gray) 

22 (note only one line for the prescribed value for CTL). 

23 Figure 9.: Zonal average TOA OLR and TOA ASW in Wm’ 2 for DJF (top), JJA (middle), and 

24 ANN (bottom). The line colors used for MAC, CTL, OBS are shown in the middle panel. 

25 Figure 10.: Distribution of simulated OLR minus OBS CRE TOA in Wm’ 2 . The rows have 

26 DJF (top), JJA(middle) and ANN (bottom). Right (left) columns are for CTL (MAC) data. 

27 Figure 11.: Same as Figure 10, but for SWA CRE TOA in Wm 2 . 

28 Figure 12.: DJF simulation for a year with MAC GCM using 8x the sea salt particles. Biases 

29 between the simulated minus observed SW CRE TOA are virtally gone. 
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2 

1 Figure 13.: Simulated ten year mean TOA net radiation (Wm'~) for DJF(top) , JJA (middle), 

2 and Annual mean (bottom) by MAC and CTL runs (left 2 panels) and MAC minus OBS and 

3 CTL minus OBS (right two panels) ; CERES data represents OBS. 
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CLOUD PHYSICS & MICROPHYSICS 


3 Figure 1 


4 Figure 1.: Cloud Physics and Microphysics of baseline GEOS-5 GCM and its modification to 

5 allow McRAS-AC to replace the relevant moist physics module of the GCM 
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Precipitation [mm/day] 

a) MAC, DJF b) CTL, DJF c) MAC-OBS, DJF d) C7L-OBS, DJF 



i) MAC, ANN j) CTL, ANN k) MAC-OBS, ANN I) CTL-OBS, ANN 



2 Figure 2 

3 Figure 2.: Simulated ten year mean precipitation (mm day' 1 ) for DJF(top) , JJA (middle), and 

4 Annual mean (bottom) in MAC and CTL runs (left 2 panels) and MAC minus OBS and CTL 

5 minus OBS (right two panels) ; GPCP data represents OBS. 
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3 Figure 3 


4 Figure 3.: Statistically significant precipitation differences (mm day-1) are color coded: MAC 

5 minus CTL for DJF (top) , JJA (middle), and Annual mean (bottom) using a 2-tailed student t- 

6 test. 
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Cloud Fraction [%] 



2 Figure 4 

3 Figure 4.: Ten year mean zonal average cloud Fractions for a) entire column atmosphere, b) 

4 400hPa to top, c) 400-700 hPa, d) and 700hPA to surface in MAC and CTL simulations 

5 versus ISSCP and MODIS Climatology. 
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1 

2 Figure 5 

2 

3 Figure 5.: Total Cloud Water Path for the full depth of the atmosphere in gm'" for DJF(top) , 

4 JJA (middle), and Annual mean (bottom) in MAC and CTL runs and MAC minus OBS (left 

5 two panels) and CTL minus OBS (right two panels); MODIS data implies OBS. 
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a) MAC, DJF b) CTL, DJF c) MAC-OBS, DJF d) CTL-OBS, DJF 
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2 Figure 6 

3 Figure 6.: Same as 3 except for Cloud Liquid Water Path; SSMI data represents OBS 
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Total Water Path [g/m 2 ] 
GOCART— OBS, JJA 



2 Figure 7 

2 

3 Figure 7.: Total Water Path for JJA in g m' ; interactive aerosol chemistry minus GOCART. 
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Effective Radius [Aim] 



2 Figure 8 

3 Figure 8. Zonal Average plots of a) ice and b) liquid cloud particle effective radius (pm). 

4 Different colors are used for Land, Ocean, and Total (colors legends are displayed in the 

5 Plots) for MODIS data (dotted), MAC simulations data (solid), and CTL simulations (gray) 

6 (note only one line for the prescribed value for CTL). 
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2 Figure 9 

3 Figure 9.: Zonal average TOA OLR and TOA ASW in Win' 2 3 4 5 for DJF (top), JJA (middle), and 

4 ANN (bottom). The line colors used for MAC, CTL, OBS are shown in the middle panel. 
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3 Figure 10 

4 Figure 10.: Distribution of simulated OLR minus OBS CRE TOA in Wm' 2 . The rows have 

5 DJF (top), JJA(middle) and ANN (bottom). Right (left) columns are for CTL (MAC) data. 
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2 Figure 11 

3 Figure 11.: Same as Figure 10, but for SWA CRE TOA in Wm' 2 . 

4 


50 


SW CRE TOA [W/m 2 ] 

SSx8— OBS, DJF 



60E 120E 180 120W 60W 


2 Figure 12 
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Figure 13 
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